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X-ray spectra of black hole binaries in the standard high/ soft state were stud- 
ies \ ied comprehensively by using ASCA GIS data, and partially RXTE PCA data. 

A mathematical disk model was applied to several black hole binaries to see if 
the observed accretion disk temperature profile was consistent with that expected 
from the standard accretion disk model. This model is called p-free disk, and as- 
sumes that the spectrum is composed of multi-temperature blackbody emission of 
which the local temperature T(r) at radius r is given by T(r) oc r~ v with p being 
a positive free parameter. The standard disk, where gravitational energy of the 
^ ; accreted matter is effectively released as blackbody radiation, roughly requires 

p ~ 3/4, while a small deviation is expected depending on the inner boundary 
conditions, general relativistic effects and disk vertical structures. Our sample 
objects included LMC X-l, LMC X-3, XTE J2012 + 381, and GRO J1655 - 40. 
During the ASCA observations, these black hole binaries showed characteristics 
of the standard high/soft state. Under the standard modeling of high-state black 
hole binaries, the sources show disk inner temperatures of 0.76-1.17 keV, the 
disk fraction to the total 0.7-10 keV flux of 54-98%, and absorption columns of 
(0.7-12) xlO 21 cm -2 . The best-fit values of p were found in the range of 0.6-0.8, 
and the standard value of p = 3/4 was accepted for all the sources. The obtained 
values of p are also compared with those expected for the standard accretion disk 
in the Schwarzschild metric by using the so-called GRAD model. The observed p- 
values were indeed found to be consistent with those expected from the standard 
accretion disk in the Schwarzschid metric. 
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Subject headings: Black hole — X-Rays: binaries — Accretion disks 

1. Introduction 

From the classical perspective, spectral states of black hole binaries are classified into 
two classes (e.g., Tanaka 1997 and references therein). One is the low/hard state in which X- 
ray spectra are represented by a single power-law with an exponential cut-off around several 
tens of keV to 100 keV. This low state is realized when the mass accretion rate is relatively 
low. The other is the high/soft state, which is known as a high accretion-rate state. In 
the high/soft state, X-ray spectra consist of a dominant soft component and a power-law 
tail. The dominant soft component is believed to be from the radiatively efficient, optically 
thick, geometrically thin disk (the "standard disk", Shakura & Sunyaev 1973) reaching the 
innermost stable orbit around a central black hole. In fact, when the soft component is fitted 
with a "multi-color disk" model (MCD model or diskbb in xspec 1 ; Mitsuda et al. 1984) 
that approximates the standard accretion disk by ignoring its inner boundary condition, the 
obtained inner radius is kept remarkably constant against any large intensity variation (e.g., 
Makishima et al. 1986; Tanaka 1989; Takizawa 1991; Ebisawa et al. 1993, 1994). Moreover, in 
many black hole binaries, including Cyg X-l (Dotani et al. 1997), once the true inner radius, 
Ri n , is estimated by the MCD fit, it is found to be consistent with R ms = 6R g (= 6GM/c 2 , 
with black hole mass M), which is the last stable orbit for a non-spinning black hole, Here, 
from an apparent inner radius, r- m , one of two fit parameters of the MCD, the true inner 
radius, i?i n , can be calculated as R{ n = r in • k 2 • £ with an inner boundary correction factor, 
£ = 0.412 (Kubota et al. 1998), and a correction factor of spectral hardening, k = 1.7-2.0 
(Shimura & Takahara 1995). 

After this standard view was established, recent spectral studies of black hole bina- 
ries have revealed that there seems to be three sub-statuses in the high/soft states: the 
standard regime, anomalous regime (or "weak very high state"; see below), and apparently 
standard regime, in the order of increasing the mass accretion rate (Kubota, Makishima & 
Ebisawa 2001; Kubota & Makishima 2004). The standard regime corresponds to the classical 
high/soft state where the inner radius is kept constant under canonical modeling of the MCD 
plus power-law. The anomalous regime is characterized by thermal inverse Compton scatter- 
ing by high-temperature electrons around the optically thick disk that extends down to -R ms , 
and by the existence of type-B QPOs (e.g., van der Klis 2004; Remillard et al. 2002). This 
anomalous regime is basically the same as a very high state (VHS; Miyamoto et al. 1991), 
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which was redefined as a steep power-law state by McClintock & Remillard (2003), and thus 
the anomalous regime was renamed as a weak VHS by Kubota & Done (2004). Though the 
soft disk component no longer dominates the hard spectral component, the standard disk 
picture is still justified for the underlying optically thick disk in this anomalous regime (Kub- 
ota et al. 2001; Kubota & Makishima 2004). Together with these three stable sub-states in 
the high state, there seems to be an intermediate case, called an intermediate state in the 
literature, or a strong VHS by Kubota & Done (2004), at an occasion of the low-to-high 
transition. 

The spectral shape of the apparently standard regime, the most luminous regime, re- 
sembles that of the classical high/soft state, and consists of the dominant soft component 
and a very weak power-law tail. However, the apparent disk inner radius, r in , obtained 
under the MCD fit is not constant, but varies as oc T~ l , where T in is the highest color 
temperature of the disk, while it is kept constant in the standard regime. Plainly, the data 
in the apparently standard regime are not consistent with a constant inner radius, constant 
emissivity and constant colour temperature correction, which are justified in the standard 
regime. We here define this inconstancy of ri n in the apparently standard regime as a de- 
viation from the simple standard picture. This deviation is now commonly seen in black 
hole binaries in the most luminous case, including GRO J1655 — 40 (Kubota, Makishima, & 
Ebisawa 2001) and 4U 1630 - 47 (Abe et al. 2005). Beyond our Galaxy and Large Magel- 
lanic Cloud (LMC), the apparently standard regime has also been suggested in off-nuclear 
ultra-luminous (3 x 10 39 ~ 10 40 erg s _1 ) compact X-ray sources in nearby galaxies (hereafter 
ULXs; Makishima et al. 2000; Fabbiano 1989). That is to say, the ULXs often show the 
MCD- type spectra, and the values of r in are not kept constant, but proportional to T m l 
(Mizuno, Kubota & Makishima 2001), like black hole binaries in the apparently standard 
regime. 

This deviation can be interpreted as a break of the standard disk (Kubota & Makishima 
2004), or can also be caused by a different color temperature correction (Gierlihski & Done 
2004). As is done for XTE J 1550 - 564 by Kubota & Makishima (2004), the break of the 
standard disk is considered in detail by investigating the local temperature gradient of the 
disk with a so-called p-free disk model (a local model for XSPEC; Kubota & Makishima 2004), 
with which the spectrum consists of the sum of blackbody emission with a local temperature 
of T(r) = T in ■ (r/r in )~ p (r is a distance from a central object). In the Newtonian metric, the 
local effective temperature of the standard disk is given as 




(1) 



with the Stefan-Boltzmann constant, a, and the mass accretion rate, M (Shakura & Sunyaev 
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1973). Figure 1 exemplifies the temperature profile of the standard disk. Thus, the p- value 
is expected to be 3/4 by ignoring the boundary condition at r m R [n . This is the case of the 
MCD model (Mitsuda et al. 1984). With another cooling process, the p- value is predicted 
to become smaller than 3/4 (e.g., Watarai et al. 2000). Figure 2 exemplifies the 0.7-10 keV 
spectra of the p-free disk model with p=0.5, 0.75(=3/4) and 1.0. This figure shows that the 
p-free disk with a smaller value of p reveals an X-ray spectrum in which the soft energy band 
is more enhanced than the standard disk of p = 3/4. 

However, as already pointed out by Kubota & Makishima (2004), even in the standard 
accretion disk, the observed value of p is expected to be slightly different from 3/4 due to 
(1) the disk inner boundary condition (see Fig. 1), (2) general relativistic effects, and (3) the 
disk vertical structure including electron scattering. Therefore, to discuss the break of the 
standard disk, knowledge of the temperature profile of the standard disk must be required. 
It is thus important to obtain an observationally acceptable range of p in well-known black 
hole binaries in the standard high/soft state (standard regime), where the radiatively efficient 
standard disk is justified. With this perspective, this paper presents a detailed study of the 
spectral shape of the standard-state Galactic/Magellanic black hole binaries observed with 
the Gas Imaging Spectrometer (GIS; Ohashi et al. 1996; Makishima et al. 1996) onboard 
ASCA. The ASCA GIS has a fine spectral capability in the range of 0.7-10 keV, and is 
free from pile-up, even for the bright (X-ray flux of ~ 1 Crab) Galactic objects. Therefore, 
the GIS is the best detector to examine the structure of the optically thick accretion disk 
around bright Galactic/Magellanic black hole binaries in detail. In this paper, we discuss 
the spectral shape by using the p-free disk model, following previous reports. 

Before the present work, information of our sample is summarized in § 2. In § 3 and § 4, 
the observed data are analyzed, and the results are summarized. In § 5, the observational 
result is briefly discussed in comparison with the standard disk picture and the latest models 
for the radiatively efficient, geometrically thin disks. 

2. Our sample and observations 

Many black hole binaries were observed with ASCA in their high/soft state (e.g., Cyg X- 
1, GRS 1009-45, LMC X-l, LMC X-3, GRO J1655-40, GRS 1915 + 105, XTE J2012 + 381, 
XTE J1550 - 564, and 4U 1630 - 47). In this paper, the ASCA GIS spectra of LMC X- 
1, LMC X-3, XTE J2012 + 381 and GRO J1655 - 40 are reanalyzed based on the p-free 
disk model. They all show characteristics of the standard high/soft state, and show an 
inner temperature of Tj n > 0.7 keV, ensuring that their Wien shoulders of local diluted 
blackbody radiation at the inner edge of the disk appear at > 2 keV. In addition, the 
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absorption columns of iV H are not too high (Nu < 10 22 cm -2 ) to study the disk spectra in 
detail. The ASCA data of XTE J1550 — 564 are not included, because they characterize 
the VHS (Kubota & Done 2004) rather than the classical high/soft state. In Table 1, the 
observation logs of these four sources are shown with the estimated central object mass, 
inclination angles, i, and distances, d. The large mass of the central objects in LMC X-3, 
LMC X-l, and GRO J1655 — 40 fiducially measured from binary motion, makes these three 
objects confirmed black holes (e.g., McClintock and Remillard 2003), and X-ray data of 
XTE J2012 + 381 characterize its nature as an accreting black hole in the classical high/soft 
state. 

For each source, the GIS events were extracted from a circular region of 6' radius centered 
on the image peak, after selecting good time intervals in a standard procedure. Average 
energy spectra were constructed from the two GIS detectors. For the two Galactic sources, 
XTE J2012 + 381 and GRO J1655 — 40, the high bit-rate data were analyzed after correcting 
the dead time in reference to the count-rate monitor data (Makishima et al. 1996). The 
medium bit-rate data were not used, because even the monitor scaler was overflown and a 
dead-time correction was not possible. A systematic error of 1% is included to each energy 
bin of the spectra, in order to account for calibration uncertainties. Because the two LMC 
sources are much dimmer, both the high bit-rate data and the medium bit-rate data are 
utilized, and the systematic error is not included. A background was completely negligible 
compared to the fluxes of the Galactic sources, and thus was not subtracted. For LMC 
sources, the background is not negligible, and is constructed from a source-free region of the 
GIS field of view. 

There were simultaneous RXTE observations corresponding to the ASCA observations 
of XTE J2012 + 381 and GRO J1655-40, on 1998 May 29 and 1997 February 26 respectively. 
The RXTEPCA data are very useful to estimate the power-law contribution, which is critical 
to examine the spectral shape of the dominant soft component with the p-free disk model. 
We used heasoft 5.2 in which systematic difference in calibration of GIS and PCA pointed 
e.g., by Done, Madejski, Zycki (2000) was corrected. Thus, for these two sources, the RXTE 
PCA data were simultaneously analyzed with the ASCA GIS data. For the PCA data 
reduction, a top layer was used from all five units, using standard exclusion criteria (target 
elevation less than 10° above the limb of the Earth; pointing direction more than 1.2 arcmin 
from the target; data acquired within 30 minutes after spacecraft passage through South 
Atlantic Anomaly). The background was estimated for each observation using the software 
package pcabackest (version 2.1e), supplied by the RXTE Guest Observer's Facility at 
NASA/GSFC. A response matrix was made for each observation by utilizing the software 
package PCARSP version 8.0. In order to take into account the calibration uncertainties, a 
systematic error of 1% was added to each bin of the PCA spectra. The 3-20 keV data were 



- 6- 



used in the present work, because there is sometimes a residual structure in the 20-35 keV 
range associated with the Xe-K edge at 30 keV. 

3. Analyses and Results 

3.1. XTEJ2012 + 381 

Responding to the discovery by the All Sky Monitor (Levine et al. 1996) on board 
RXTE, ASCA observed XTE J2012+381 on 1998 May 29-30 with a net exposure of 38 ks 
(White et al. 1998). As previously shown by White et al. (1998), the energy spectrum was 
very soft, and was successfully described with the canonical MCD plus power-law model. 
Together with the high flux and low fractional contribution from the power-law tail, this 
fact makes the object ideal for the present study. The RXTE pointing observations of this 
source were performed 24 times; Vasiliev, Trudolyubov & Revnivtsev (2000) reported that 
the inner disk radius is kept constant under the canonical MCD plus power-law model. One 
of the RXTE observations on 1998 May 29 17:48:00-18:46:56 (UT) was simultaneous with 
the ASCA observation. The GIS spectrum was reconstructed by using the data during this 
window with an exposure of 2.4 ks. The 3-20 keV RXTE PCA spectrum and the 0.7-10 
keV ASCA GIS spectrum were jointly analyzed. 

The overall spectrum was first examined with the MCD plus power-law model. Except 
for a normalization factor, all of the model parameters were constrained to be the same 
between the GIS and the PCA data. Both these components were assumed to undergo 
absorption with a common column density, Ah- The model reproduced the simultaneous 
data very well with x 2 / u = 152.1/152 on the condition that a Fe-K edge around 7.6 keV 
was included in the power-law component. The Fe-K edge structure was significant in terms 
of the F-test (F(3, 152) = 218), and the smeared-edge model (smedge in XSPEC; Ebisawa 
et al. 1994) gave an acceptable fit to the structure. Figure 3a shows the raw data and the 
residuals between the data and the model. In Table 2, the best-fit parameters are shown 
together with the fraction of the disk luminosity to the total luminosity in the range of 
0.7-10 keV. Though the values of d and i are not known, the true inner radius, R in) and 
the disk bolometric luminosity, Ldisk? were calculated with the assumptions of d — 10 kpc 
and i = 60°. Here and hereafter, a spectral hardening factor of k = 1.7 and a boundary 
correction factor of £ = 0.412 are used to estimate of R in . 

The successful fitting of the canonical model means that the temperature gradient does 
not conflict with p = 3/4, which is the case of the MCD. However, to obtain both the best-fit 
value and an acceptable range of p, the soft component was next represented with the p-free 



-7- 



disk model instead of the MCD model (hereafter model A). Other fitting conditions were 
the same as in the canonical fitting. As also shown in Table 2, the best-fit p was determined 
as 0.71 ± 0.02 within 90% confidence. Even though the power-law parameters have been 
accurately constrained thanks to the RXTE data, our basic assumption that the hard com- 
ponent has a single power-law shape down to ~ 1 keV would not be warranted; the value of 
p might change if these systematic uncertainties were taken into account. This uncertainty 
was considered by applying an additional absorption colmun, iV^ ard , to the power-law com- 
ponent alone (hereafter model B). The results are also summerized in Table 2. The value of 
p is obtained as 0.8 ± 0.1 with large iV£ ard of 2.5 ± 0.8 x 10 24 cm" 2 . The bottom two panels 
in Fig. 3a show the residuals between the data and models A and B. 

Figure 4 shows unfolded spectra with the best-fit predictions of models A and B. This 
figure shows that the power-law contribution becomes maximum with model A, while it 
becomes minimum with model B, and thus the fittings with these two models give the most 
conservative range for p. The acceptable range of p is thus determined to be 0.7-0.9, even 
when taking into acount the systematic uncertainties of the hard component. Hereafter, all 
of the sources are considered based on with the canonical model, model A and model B. 



3.2. GROJ1655 -40 

The X-ray transient GRO J1655 — 40 (Nova Sco 1994) is a secure black hole binary 
system with dynamical mass evidence of M — 6.0-6.6 M that exhibits superluminal jets 
(e.g., Hjellming & Rupen 1995; van der Hooft et al. 1998; Shahbaz et al. 1999; Hannikainen 
et al. 2000). Based on multiple RXTE pointings of the 1996-1997 outburst of this source, 
Kubota et al. (2001) showed that it exhibited the previously described three spectral regimes 
of a high luminosity state with the critical disk luminosity of 1-2 xlO 38 erg s _1 . 

ASCA observations of GRO J1655— 40 were performed five times from 1994 to 1997. The 
source was found in the high state in 1995 August 15-16 and 1997 February 25-28 (Zhang et 
al. 1997; Ueda et al. 1998; Yamaoka et al. 2001). Referring to the spectral parameters of the 
1996-1997 outburst (Kubota et al. 2001), on one hand, the source was found on the highest 
end of the standard regime during the 1997 observation. On the other hand, the spectral 
parameters of the 1995 observation (e.g., Zhang et al. 1997) reveal the characteristics of the 
anomalous regime (weak VHS), and thus the strong inverse Comton scattering may change 
the spectral shape from direct disk emission. Thus, we focus our discussion on the 1997 
observation. The ASCA observation of 1997 February 26 was simultaneously performed 
with RXTE (e.g., Yamaoka et al. 2001). Net exposures of 2.4 ks of the ASCA GIS high- 
bit rate data and 5.4 ks of the RXTE PCA data were obtained after the standard data 
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reductions. 

Following the case of XTE J2012 + 381, the 0.7-10 keV GIS spectrum and the 3-20 
keV PCA spectrum were simultaneously analyzed with the canonical MCD plus power-law 
model, model A, and model B. All of the models assume that the soft disk component (MCD 
or p-free disk) and the power-law components are commonly modified by the inter-stellar 
absorption and other absorption features. The fit results are shown in Table 3 and Fig. 36. 
In this spectral fitting, a negative Gaussian and three-edge structures were included. By 
referring to Yamaoka et al. (2001), the center energy of the negative Gaussian was fixed at 
6.81 keV, and the energies of the first two edges were fixed at 7.71 keV and 8.81 keV. The 
obtained values of the equivalent width of the negative Gaussian, and the optical depth of 
the edges are consistent with those shown in Yamaoka et al. (2001). In our joint fit, the 
other edge structure was required at ~ 10.4 keV of great significance (F(2, 150) = 156.9 in 
the case of the canonical model). Due to the high temperature of the disk, and very low 
intensity of the power-law component (see Fig. 36) , the value of T of the power-law could not 
be precisely determined even with the PCA data. Therefore, under a fitting with models A 
and B, the value of T was limited to the range of 2.0-2.5, which is the typical value for the 
power-law component in the standard high/soft state. 

With the canonical fit, the true inner radius was obtained as i?i n = 22. 6^03 km, for 
d = 3.2 kpc and i = 70°. The disk luminosity was obtained as Ldisk = 9 x 10 37 erg s _1 , which 
is roughly 10% of the Eddington luminosity of a 6.0-6.6 M black hole. The confidence 
ranges of i? in and Ldisk given in Table 3 include systematic uncertainties of the estimation 
of d together with statistic errors. With models A and B, the acceptable range of p was 
obtained as 0.68-0.73. 



3.3. LMC X-l 

One of the first identified X-ray sources in the Large Magellanic Cloud, LMC X-l 
(Johnston, Bradt, & Doxsey 1979), is known as a secure black hole binary with M=4-10 M & 
and i < 63°. This object was observed with ASCA on 1995 April 2-3 with a net exposure of 
21.7 ks. Makishima et al. (2000) reported that the ASCA GIS spectrum was well reproduced 
by the canonical MCD plus power-law model. The same data were reanalyzed based on the 
same condition for the reduction in Makishima et al. (2000). The fitting results with the 
canonical model, model A, and model B are given in Table 4 and in Fig. 3c. Same as the 
case of XTE J2012 + 381, the smedge is significant with F(3, 110) = 31.1. With d = 50 kpc 
(Freedman et al. 2001) the true inner radius was obtained as i?i n = 59±2- (cosi/ cos30°) -1 / 2 
km from the canonical model. For a 4-10 M Q black hole, this is consistent with ~ 6R g . The 
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calculated value of L disk , ~ 1.5 x 10 38 • (cos?/ cos30°) _1 erg s" 1 , is 9-20% of the Eddington 
limit for a 4-10 M Q black hole with an inclination of % < 63°. The acceptable range of p 
was determined to be 0.7-1.1 with models A and B. 



3.4. LMC X-3 

Together with LMC X-l, LMC X-3 is a persistent black hole binary in LMC with 
M=5.9-9.2 M . It was observed with ASCA twice on 1993 September 22-23 and 1995 
April 14-15, with net exposures of 20.6 ks and 18.4 ks, respectively. The ASCA results 
have already been reported by Makishima et al. (2000). The same data as in Makishima 
et al. (2000) were reanalyzed by applying the p-free disk model in this paper. The best-fit 
parameters based on the canonical model are given in Table 5. The raw data are shown with 
the best-fit canonical model in Fig. 3d-e. Without any edge structure, the canonical model 
has given acceptable fits to 0.7-9 keV GIS spectra, with ^jv = 102.2/109 and 103.8/109, 
for 1993 and 1995 observations, respectively, and the fitting parameters are consistent with 
those from Makishima et al. (2000). With d = 50 kpc and i = 66°, the innermost radii and 
the disk bolometric luminosities were calculated to be Ri n = 47 ± 2 km and 46 ± 2 km, and 
Ldisk = 8.5 x 10 37 erg s _1 and 1.7 x 10 38 erg s -1 for each observation. The obtained R[ n are 
consistent with 6R g , and the Ldisk values are roughly 6-10% and 12-19% of the Eddington 
luminosity of 5.9-9.2 M Q BH. The p-free disk model was also applied to this source as 
model A and model B, and the acceptable ranges of p were determined to be 0.51-0.75 and 
0.60-0.70, for the 1993 and 1995 observations, respectively. 

4. Discussion 

4.1. Summary of the observational results 

The disk temperature profiles of four black hole binaries (XTE J2012+381, GRO J1655- 
40, LMC X-l and LMC X-3) were examined with the p-free disk model. The uncertainties 
of the contribution of the hard power-law component were considered by fitting the data 
with two cases, model A and B. In Fig. 5, the best-fit values of p under model B are plotted 
against those under model A for our sample. The obtained values of p are almost the 
same between these two cases. Hereafter, an average of the best-fit values of p with model 
A and B is applied as the most reliable value with its error, which covers an acceptable 
range of p with these two models. In Fig. 6a, the resultant p-values are plotted against 
Ti n . Here, T in under the MCD fit is applied instead under the p-free disk fit, to escape 
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from any systematic coupling between the parameters of the p-free disk model. For MCD 
temperatures of T in =0.76-1.17 keV, the central values of p are found to be between 0.6 and 
0.8 with an average of 0.70 ± 0.05; the standard p- value of 3/4 is accepted for all sources. 

Figure 6b-c show other scatter plots, where the obtained values of p are plotted against 
the 0.7-10 keV disk fraction and Nu- The same as in Fig. 6a, the values of the disk fraction 
and iVn are based on the canonical MCD plus power-law model. In this figure, there is no 
dependence of p on the other spectral parameters. It is thus concluded that the observed 
0.7-10 keV GIS spectra of the standard high/soft state black hole binaries are actually very 
consistent with the simple standard accretion disks of p = 3/4 for T[ n =0.76-1.17 keV, the 
disk fraction of 54-98%, and A" H =0.7-12x 10 21 cm" 2 . 

To check the p-dependence on the disk luminosity, Eddington ratio, and black hole spin, 
in Fig. 7 the values of p are plotted against i, Ldisk, -^disk/^Edd) and R in /R g . The result 
of XTE J2012 + 381 is not shown in this figure because of a lack of optical information. 
Together with 90% statistical errors, the systematic uncertainties of % and d are considered 
to estimate the errors of L disk , L disk /L Edd and i? in /i? g . This figure also does not show any 
dependency of the p-values on the other parameters. 



4.2. On the inner boundary condition and the general relativistic effects 

The standard accretion disk is the basic model for a radiatively efficient, geometrically 
thin disk. A p- value of 3/4 is approximately expected for the standard accretion disk. 
However, the combination of the boundary condition of the inner edge of the disk, the 
general relativistic effects, and the disk vertical structures was expected to change the p- 
value from 3/4 near the inner edge. Therefore, it is meaningful to compare the observed 
p-values with those predicted theoretically for the standard disk. 

As described in § 1, the disk inner-boundary condition of friction free makes the tem- 
perature profile flatter than r _3//4 (i.e., the value of p becomes smaller than 3/4; see Fig. 1). 
If the Schwarzschild potential is applied instead of the Newtonian potential, the difference of 
the temperature profile around the inner boundary becomes more significant. A general rela- 
tivistic accretion disk model (grad model in the XSPEC; Ebisawa, Mitsuda & Hanawa 1993; 
Hanawa 1989) calculates the local effective temperature of the standard accretion disk in the 
Schwarzschild potential as 

3GMM f i 3i? g \ _1 f fm^ Fm^ ^(l + yfiRjr 1-1/ V / 2 XV ' 



T eff (r) = 

8vrcrr 3 V r J { V r V 4r \l-^3R~Jr l + l/y/2 

(2) 

In Fig. 1, this temperature profile is shown together with that in the Newtonian metric 
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via Eq. (1). The figure shows that the flattening of the temperature profile at the inner 
boundary is emphasized under the Schwarzschild potential. 

In addition to the difference of T eff (r), the gravitational redshift and the Doppler boost- 
ing are the main relativistic effect; the latter appears most clearly in the dependence of 
the observed spectra on the inclination angle, i. In order to test both the inner boundary 
condition and the general relativistic effects in the Schwarzschild metric, the p-values are 
calculated based on the grad spectra under the 0.7-10 keV GIS response matrices. In de- 
tail, many GRAD spectra are simulated for several conditions of % and M, and are fitted with 
the p-free disk model and the MCD model. To compare the observational results shown in 
Fig. 6a, we plot the predicted p- values for the GRAD model against T in for several conditions 
of % in Fig. Sa-b. The uncertainty of N R is ignored and considered in panel (a) and (6), 
respectively. The trends of p-dependence on T in are the same between these two panels, 
though the absolute p-values in panel (6) are distributed more widely than those in panel 
(a). In the case of face-on systems with i = 0° ~ 30°, panel (6) shows that the values of 
p are expected to appear in the range of 0.68-0.75, and do not depend on T- m very much. 
On the contrary, those for edge-on systems with % = 30° ~ 70° depend on T m significantly, 
because the gravitational boosting at the inner boundary becomes efficient. In this case, the 
values of p are distributed between 0.58-0.7. 

These results show that the relativistic effects under the Schwarzschild metric and the 
boundary condition make p smaller than 3/4 up to Ap m —0.2. By comparing Fig. 6a and 
Fig. 86, in Fig. 9, the observed values of p are plotted against those predicted for each binary 
system based on the GRAD model. Here, the GRAD predictions were obtained based on the 
observed MCD temperature and the optical information of % for each source by using Fig. 86. 
For XTE J2012 + 381 and LMC X-l without (or with poor) information on i, the errors of 
p- values for the GRAD prediction cover the range for % = — 70° and % = — 63°, respectively, 
and their central values are just taken from the center for the confidential range. 

That the observed best-fit values of p are consistent with those calculated for the GRAD 
model, and a trend toward a slightly larger value of the observed-p than the GRAD-p (Ap ps 
+0.06) is marginally suggested. 

4.3. On the disk vertical structure 

The consistency of the p-values between the observations and the GRAD calculations 
means that the disk vertical structure, which is not included in the GRAD model, should not 
change the spectral shape very much. Therefore, the assumption concerning the sum of the 
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local diluted blackbody radiations with a radially constant hardening factor, k, is justified 
for black hole binaries in the standard high/soft state. Otherwise, the vertical structure can 
shift the p-value to slightly larger, as marginally suggested in Fig. 8. 

Is this observational result consistent with the current models for a geometrically thin 
disk that includes the disk vertical structure? Even within the standard disk picture, electron 
scattering in the disk direction is expected to be important. Since electron scattering under 
optically thick conditions usually becomes more important for higher temperature (e.g., 
Rybicki & Lightman 1979), emission from the higher temperature inner part of the disk may 
be suppressed more than the lower temperature outer part. Thus, with electron scattering, 
the p- value roughly becomes smaller than the standard value of 3/4. The real disk is much 
more complicated, and thus such a brief estimation may not be realistic. In fact, the disk 
spectra that take account for electron scattering have been computed by many authors, 
including Shimura & Takahara 1995, Ross & Fabbian (1996), and Merloni, Fabbian & Ross 
(2000), on many conditions. 

Recently, significant theoretical advances have been made in understanding the nature 
and the vertical structure of the optically-thick and geometrically-thin disk. Turner (2004) 
shows that magnetic stress should play an important role in the disk, and is considered to 
make the radiatively efficient disk stable. Bound-free metal opacity is also thought to be 
important (Davis et al. 2005). Together with these theoretical studies, the observationally 
determined temperature profile for the disk can be a good probe to see the distortion of the 
spectral shape due to the disk vertical structure. 

We are grateful to T. Shimura for comments on the disk model by Shimura & Takahara 
(1995), and A. K would like to thank C. Done, S. Mineshige and Y. Terashima for valuable 
discussions. A.K. is supported by a special postdoctoral researchers program in RIKEN. 
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Table 1. Observations with ASCA and optical data of our sample. 



object 


BH mass 


inclination 


Distance 


observation date 


exp. 


reP 




Mq 


[degree] 


[kpc] 


yyyy/mm/dd 


[ks] 




XTE J2012 + 381 








1998/5/29-30 


2.8 


1,2,3 


GRO J1655 - 40 


6.0-6.6 


69.5 ± 0.08 


3.2 ±0.2 


1997/2/26 


2.4 


4,5,6,7,8,9,10 


LMC X-l 


4-10 


< 63 


55 


1995/4/2-3 


21.7 


11,12,13,14 


LMC X-3 


5.9-9.2 


64-68 


55 


1993/9/22-23 


20.6 


11,14,15 










1995/4/14-15 


19.9 





a References — 1 White ct al. 1998; 2 Vasiliev et al. 2000; 3 Hynes et al. 1999; 4 Hjellmiung & Ruppen 
(1995); 5 Orosz & Bailyn (1997); 6 Shahbaz ct al.(1999); 7 Zhang ct al. (1997); 8 Ueda et al. (1998); 
9 Yamaoka et al. (2001); 10 Kubota ct al. (2001) n Fcedman ct al. (2001); 12 Hutchings et al. (1987); 
13 Cowlcy ct al. (1995); 14 Makishima et al. (2000); 15 Cowley et al. (1983); 



Table 2: Best-fit parameters obtained by a simultaneous A SCA/ RXTE observation of XTE 
J2012 + 381. 



parameters 


MCD 


p-free 


disk 






A 


B 


ZV H (10 22 cm- 2 ) 


1 24+ - 04 


1.19 ±0.05 


1.1 ±0.1 


MCD or p-free disk 








fcHn (keV) 


0.761 ±0.005 


0.77 ±0.02 


0.76 ±0.01 


norm. (xl0 3 ) a 


1.62 ±0.05 


1 4+ - 4 


2.ol8:i 


P 


(0.75) 


0.71 ±0.06 


0.8 ±0.1 


power-law 








^yhard ( 1Q 22 cm -2) 






25 ±8 


r 


9 o R +0.06 


2.3 ±0.1 


2.6 ±0.1 


norm. 6 


0.44 ±0.07 


0.4 ±0.1 




Absorption edge 








E (keV) 


7 6+ - 1 


7.3 ±0.3 


o O+0.2 
X-3-0.3 


T 


> 0.9 


2+_\ 


> 0.8 


width (keV) 




c + 19 
°-2 


9 fi +0.4 




152.1/152 


151.1/151 


139.7/150 


disk fraction c 


0.66 


0.67 


0.95 


Rin [km] de 


68 ± 1 






L disk [10 38 erg s^] d 


1.4 







Note — Errors represent 90 % confidence limits. 
a In unit of (r; n /km) 2 • cosi ■ (d/10 kpc) -2 
b In unit of photons s -1 cm -2 keV -1 at 1 keV 

c A ratio of flux of the MCD component to the total flux in 0.7-10 keV. 
d Estimated for assumed distance of 10 kpc and inclination of 60°. 
e With k = 1.7 and £ = 0.41 
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Table 3: Same as Table 2, but for GRO J1655-40. 



parameters 


MCD 


p-free 


disk 






A 


B 


iV H (10 21 cm" 2 ) 




7 1+0.2 
0.3 


7.1 ±0.3 


disk 








kT in (kcV) 


i 171 +0.005 

i - 1 ' s --v.om 


1.19 ±0.01 


1 OQ+0.01 

l.ZU_Q Q 2 


norm. (xlO 3 ) 


i 91+0.02 

-"^ -0.04 


1.0 ±0.1 


97+ - 07 

u - a ' -0.13 


P 


(0.75) 


0.71±°;«I 


70+ 03 


power-law 








^hard (- 10 24 cm -2) 






7+8 


r 


9 7+0.3 

Z -'-0.5 


(2.2>t° ; 3 2) 


(2.5>+8. B ) 


norm 


o.5t 3 ; 2 


n+ 018 


S5+ 07 


Absorption structures" 








edgel h t 


0.15 ±0.03 


0.19 ±0.04 


o.is±8:S 


edge2 h t 


u - zi -0.05 


n 99+O.O4 
u ' zz -0.05 


0.21^ 


edge3 E [keV] 


10.4 ±0.2 


10.4t° i 


10-4±C;1 


T 


0.27 ±0.05 


o.29ig;8f 


n 90+0.05 

U - z8 -0.08 


gaussian c EW [eV] 


61 ±18 


70+16 

' °-23 


79 ± 18 


X 2 /v 


123.1/150 


116.4/149 


114.8/148 


disk fraction 


0.975 


0.990 


0.999 


R m [km] de 


22.6 ±0.2 t\i 






L disk [10 38 erg s" 1 ] # 


0.88 ±0.11 






a Referring to Yamaoka < 


3t al. (2001). 






b Energies of edgel and edge2 are fixed at 


7.71 keV and ; 


3.81 keV, 



respectively. 

c Negative gaussian to represent an absoroption line. Line center 

energy and a are fixed at 6.81 keV and 160 eV, respectively. 
d Central values are estimated for d = 3.2 kpc and i = 70°. 
e Errors include the statistic error (first one) and the uncertainties of 

i and d (second one). 
f Errors include the uncertainties of i and d. 
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Table 4: Same as Table 3, but for an ASCA observation of LMC X-l. 



parameters 


MCD 


p-frcc disk 








Q 


N R (10 21 cm" 2 ) 


6.0+ 2 , 

— o 


a n+ll.2 
b - U -0.5 


4 - 4 -0.4 


disk 








fcT in (keV) 
norm 


0.830™ 

— U.UUO 

85+^ 




w ^ — 0.02 

103^ 4 


0.82 ± 0.04 
112+1? 


P 


(0.75) 


0.8 ±0.1 


0.8+°? 

—0.1 


power-law 








^hard ( 1Q 22 ^-2) 

r 

norm 


— 
2.4 ±0.2 
1 1 + - 04 

U.ll_o.02 


(2.5>±g. B ) 

n U++ - 1 

u - 14=c -0.07 


0-5+05 
(2-5>±8. 6 ) 
0-1371X 


Absorption edge 








E (keV) 


7 O+0.4 


7.1 ±0.4 


7.1 ±0.4 


T 


> 0.6 


> 0.5 


> 0.5 


width (keV) 


2± 1 


2±1 


2± 1 


r> 


94.6/110 


93.7/109 


92.8/108 


disk fraction 


0.68 


0.64 


0.69 


Rin [km] a 


59 ± 2 ±f 






i disk [10 38 erg s" 1 ]* 


1 5+ 1 - 4 






a Central values are estimated for d 


= 50 kpc and 


i = 30°. 



Table 5: Same as Table 3, but for ASCA observations of LMC X-3. 







1993 






1995 




parameters 


MCD 


p-free disk 


MCD 


p-free disk 






A 


B 




A 


B 


iV H (10 21 cm" 2 ) 


1.0 ±0.4 




9+ ul 


7+ - 4 

u -'-0.5 


4+ ' 1 


3± 2 


disk 














kT in (keV) 
norm 


0.81 ±0.2 
25 ±3 


o-83lS:g 

19 ±8 


0.83 ±0.05 

W +15 


°- 97 -o.i 
24 ±2 


98+ 04 

u ' yo -0.02 

21 ±6 


0.99 ±0.04 

let! 


P 


(0.75) 


0.671H 8 


u - oo -0.04 


(0.75) 


O.Tllfcg? 


u - D,3 -0.03 


power-law 

^yhard (- 10 22 cm -2) 

r 

norm. [xl0~ 2 ] 


9 7+0.1 

6±2 


(2-5>±8. 5 ) 
4±S 


3+? 
(2.5> + ^) 


9 7+O.2 
z -'-0.3 

5±2 


(2-5>±g. B ) 

o+l 

J -2 


3+? 
(2-5>±g. B ) 

q+6 

J -2 




102.2/109 


102.6/108 


98.9/107 


103.8/109 


104.0/108 


101.4/107 


disk fraction 


0.54 


0.63 


0.86 


0.76 


0.82 


0.94 


Rin [km] a 


47 ±2 ±2 






46 ±2 ±2 






L disk [10 38 erg s" 1 ]* 


U - OO -0.04 






1.7 ±0.1 







a Central values are estimated for d = 50 kpc and i = 66° . 
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Fig. 1. — Radial profiles of T eff shown based on the standard accretion disk under the 
Newtonian potential (dashed line) and the Schwarzschild potential (solid line). 
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1000 g-r 

>100 r 




0.1 — ' ^ 1 1 1 1 1 1 

12 5 10 

energy (keV) 



Fig. 2. — Calculated spectra of the p-free disk model. The cases of p = 0.75 (solid line), 0.5 
(dashed line) and 1.0 (dash-dot line) are shown under the assumption of kT [n = 0.8 keV and 
(r in /km) 2 cosi/(rf/10kpc) 2 = 10 4 . The ratios to each models to the disk model of p = 0.75 
are also shown. 
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Fig. 3.— GIS and PCA spectra of XTE J2012 + 381 (a), GRO J1655-40 (6), LMC X-l (c), 
LMC X-3 in 1993 (d) and LMC X-3 in 1995 (e). The best-fit canonical MCD plus power-law 
model are superposed on the data. The residuals between the data and the canonical model, 
model A and model B are shown in the bottom three panels. 
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(a) XTE J2012+381 (model A) 



(b) XTE J2012 + 381 (model B) 




Energy [keV] 



Energy [keV] 



Fig. 4. — Unfolded spectra of XTE J2012 + 381 with the best-fitting model A (panel a) and 
model B (panel b). 



0.9 



0.8 



a) 0.7 - 
o 



O XTE J2012+381 

• GRO J1655-40 

* LMC X-l 
□ LMC X-3 



0.6 



0.5 - 



0.5 0.6 0.7 

p (model A) 



0.8 



0.9 



Fig. 5. — Best-fit values of p under model B are plotted against those under model A. 
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(a) 



(b) 



O XTE J2012+381 

• GRO J1655-40 

* LMC X-l 
□ LMC X-3 



ft 0.8 



-ft 0.8 - 



0.6 - 



0.6 - 



0.6 



0.8 



TjkeV] 



(c) 



ft 0.8 - 



0.6 - 



0.1 



1.2 



0.5 0.6 0.7 0.8 0.9 

disk fraction (0.7-10 keV) 



N H [10 21 cm 2 ] 



10 



Fig. 6. — Result of fitting with the p-free disk model. The best-fit p- values were plotted 
against kT in (a), the 0.7-10 keV disk fraction(b), and iVn(c). In panel (a), T in under the 
MCD fit is applied instead under the p-free disk fit, to escape any systematic coupling 
between the parameters of the p-free disk model. The average of the best fit values of p with 
model A and model B are plotted with marks. Each horizontal axis is based on the MCD 
fit. 
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(a) 



(b) 











• 


GRO J1655-40 


1 








* 

□ 


LMC X 
LMC X 


-1 

-3 


0.8 






























-E 




0.6 










-E 






(c) 


20 


40 

i [degree] 




60 


8C 



a O.E 



0.6 



0.1 0.2 

disk dd 



a 0.8 



0.6 





(d) 



a 0.8 



0.6 



0.5 1 1.5 

L d,sk [1° 38 ^s/s] 



10 



Fig. 7. — Same as Fig. 6, but the p-values are plotted against i (a), Ldisk (b), Ldisk/^Edd (c), 
and R in /R g (d). The errors of the horizontal axes include the systematic uncertainties of i 
and d together with the statistic errors. 
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(a)No absorption column (b)With absorption column (N H =10 21 cm 2 ) 
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• 


p=3/4 - 
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0.7 ( 
i 


> o o o o o 2 § § $ 5 8 . ft. $-f 




0.7 

i 


o o „ *?»•.., 

o O o o O O o o O O o o o o 

. . A A A A 

*a"-* A U^ 

. . • ■ • 










a 






0.6 


• i = 


□■— i=40 


0.6 


. • " • i = 


Q-- i=40 




- 1= 10 


i = 50 




■ 1= 10 


-— i = 50 




o i = 20 


■-- i=60 




o i = Z0 


■-- i=60 




A- i = 30 


i = 70 




* i = 30 


i = 70 


R 














Tm [keV] Tin [ ke V] 



Fig. 8. — Expected values of p against T in for the GRAD spectra with several inclination 
angles, i. The horizontal axis, T in , is based on the MCD fit. The uncertainty of the systematic 
coupling of the absorption column and the j9-value is ignored in panel (a) and considered in 
panel (b). The detailed criteria for the simulation are as follows. Panel(a) — (l)Many GRAD 
spectra without N H are calcurated by changing both % and M. The simulated spectra are 
convoluted with the response function of the 0.7-10 keV GIS by including 1% systematic 
errors. (2) The faked GRAD spectra were fitted with both the MCD model and the p-free 
disk model without A^h- (3) For the same GRAD spectra, the best-fit p of the p-free disk 
model are plotted against the best-fit T in of the MCD model. Panel (b) — Same as panel (a), 
but the spectra were simulated and fitted by including iVn- The plot is the result for A^h of 
10 21 cm~ 2 . This plot is not influenced by changing A^h in the range of 10 20 ~ 22 cm~ 2 . 
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0.5 0.6 0.7 0.8 

p (GRAD) 

Fig. 9. — Obtained values of p plotted against the expected values of p for each binary 
system based on the GRAD model in Fig. 8 with information of % and T ni . The errors of 
p of the GRAD expectation include the uncertainty of the estimation of %. Especially for 
XTE J2012 + 381, the range of p under the GRAD prediction is justified for % of 0-70°. The 
solid line and dashed line mean p(observed) = p(grad) and p(observed) = p(grad)+0.06, 
respectively. 



